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| Abstract. Based on a high resolution cosmological n-body simulation, we track the hierarchical 

. growth of black holes in galaxy clusters from z = 20 to z = 0. We present a census of black holes as 

function of redshift and will determine their mass assembly history under a variety of assumptions 
regarding the importance of gas accretion in black hole growth, from early supercritical Eddington 
accretion to gas-poor hierarchical assembly. Following a galaxy merger, black holes are expected 
to form, inspiral and merge after strongly radiating energy via gravitational waves. For each binary 
, black hole inspiral and merger, we determine the expected gravitational wave signal for the Laser 

Interferometer Space Antenna (LISA), and calculate the LISA event rate as a function of time. We 
will calculate the black hole mass assembly history for several black hole growth scenerios, so that 
we can explore tests to characterize each model observationally. In particular, we will study how 
well LISA observations will be able to distinguish between these very different assembly scenarios. 
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MOTIVATION 

In current theories of structure formation, the first generation of star formation seeds 
intermediate mass black holes deep within the potential of forming galaxies (Madau & 
Rees 2001, Abel et al. 2002, Heger et al. 2003). Mass of the structures that can host 
first massive black holes is limited by various feedback processes. From Population 
III supernova studies, the formation of first stars stops at z~12 (Wise & Abel 2005). 
These black holes can grow in mass by accreting gas, and as the event horizon grows 
in size, stars and stellar remnants. Later, galaxy mergers provide an impetus for these 
growing black holes to meet and form a bound system (Volonteri et al. 2003, Islam 
et al. 2003). During a galaxy merger, each black hole sinks to the center of the new 
galaxy potential due to dynamical friction, and eventually becomes bound as a binary. 
Dynamical friction then continues to shrink the orbit until the binary is hard (i.e., the 
separation between each black hole, «bbh 5 is such that the system tends to lose energy 
during stellar encounters)(Madau & Rees 2001, Holley-Bockelmann et al. 2001, Vine 
& Sigurdsson 1998). Thereafter, further decay is mediated by 3-body scattering with 
the ambient stellar background until the binary becomes so close that the orbit can lose 
energy via gravitational radiation. In studies of a static, spherical potential, it may be 
difficult for stellar encounters alone to cause the binary to transition between the 3- 
body scattering phase and the gravitational radiation regime; however in gas-rich or 
non-axi symmetric systems the binary hardens efficiently into one that emits copious 
gravitational radiation. Thereafter, it presumably coalesces. Depending on the black 
hole mass, the final stages of coalescence emit so much gravitational radiation (Thorne 




FIGURE 1. la-left- Merger tree from z = 20 to z = 3. Different halos are represented with different 
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colors. The sizes of the circles are set by the virial radius, which is related toM h ^ lo - j .. By our dry merger 
assumption, Mgni, is also related to Mh a io,viriai- Given this merger tree, we will test different models for 
black hole accretion, lb-right- 2D density projection of the simulation box with primary halo. Overploted 
are paths of dark matter halos which grow in mass (circles represent scaled virial radii) and sink toward 
primary. Only main branches of the merger tree are presented. 



1995, Cornish & Levin 2002, Vecchio et al. 2004) that they are extremely likely to be 
detected by the Laser Interferometer Space Antenna (LISA), a planned NASA mission 
to detect gravitational waves, set to launch in the next decade. Current estimates of the 
total gravitational wave signal from the cosmological growth and merger lifecycle)of 
black holes have been made using semianalytic models of merger trees, with analytic 
prescriptions for the black hole merger timescales within the analytic halos and models 
for gas accretion. In this study, we use high resolution N-body simulations to track the 
seed black holes themselves as they sink within a galaxy halo. We will then use this to 
more accurately model different black hole growth scenerios, which will yield different 
gravitational wave signatures. Ultimately, we will determine if LISA observations can 
be used as a tool to discover how black holes grow. 



GROWING SUPERMASSIVE BLACK HOLES 

In our numerical simulations we use GADGET to evolve a comoving 14.3 Mpc 3 section 
of a ACDM universe (^m=0.3, Q a =0J and h=0.7) from z = 40 to z = 0. We refine a 
sphere of 2 Mpc in the box to simulate at higher resolution with 4.9 x 10 6 high-resolution 
particles (softening length 2 kpc comoving). The rest of the box has 2.0x10 low- 
resolution particles (softening length 4 kpc comoving). The mass of each high resolution 
particle in this simulation is 8.85 x 10 5 M Q and the mass of each low-resolution particle 
is 5.66 x 1O 7 M . We use P-GroupFinder to define dark matter halos and to identify black 
holes as the most bounded particles in their host halos. 

Only dark matter halos with mass M VI >> 10 8 M Q [( 1 +z) / 10] ~ L5 can host first stars and 



FIGURE 2. 2a-Left- Number of IMBHs mergers as a function of redshift. 2b-Middle- Average ratio 
of merging black holes as a function of redshift. 2c-Right- Average of more massive components in the 
merger as a function of redshift. 



the first stars stop forming at z~12. By tracking the positions of over 3800 seed black 
holes and their host halos throughout the simulation (Fig. 1), we will model different 
black hole merger scenerios and track the growth of the black hole mass per halo as a 
function of time and environment. In particular, we will parameterize the degree of gas 
accretion involved with growing the black hole and the efficiency of binary coalescence. 

In this first attempt, however, we model the seed black hole merger timescale and 
the mass accretion history in the simplest possible way. We assume that the hardening 
timescale is rapid, that no black holes are ejected, and that there is no gas accretion 
involved. This gives us a fiducial model with which to compare more realistic black hole 
growth scenerios. 

DRY MERGER AND MERGERS PLUS ACCRETION 

We start with 200M Q initial black hole mass. Given this simple growth scenerio, we find 
a large number of mergers at each redshift, with the largest number of mergers at z ~ 6. 
At later times, 100 : 1 mergers are much more common, as Mbhi grows larger while also 
accreting smaller mass black holes from satellite halos (Fig. 2). At this stage we have 
analyzed the data in range 3<z<20. The intensity of mergers of dark matter structures 
demands more detailed analysis, to be addressed in the future work. Notice that dry 
mergers build massive black holes M>1O 6 M by redshift z=3. We add accretion to the 
previously described growth of massive black holes through mergers. For a sustained 
Eddington accretion of baryon matter the mass growth rate is Mg# / ts a i, where ts a i is 
Salpeter time-scale, ts a /~4x 10 7 yr from the recent observations. Accretion is triggered 
by mergers of dark matter halos and lasts tsai- During this time black holes doubles its 
mass. 

EXTRACTING GRAVITATIONAL WAVES 

We assume the black hole binary only emits gravitational radiation during the 'gravita- 
tional wave' phase, where the orbital separation of the binary decays most rapidly due 



FIGURE 3. 3a-left- Mass function for merging black holes. Red - redshift z=19; Green - redshift z=10; 
Blue - redshift z=3. 3b-right- Mass function for massive black holes which grew through mergers and 
accretion. There is a number of supermassive black holes at redshift z=3. 



to gravitational wave emission. In reality, this is a lower limit to the amount of gravita- 
tional radiation emitted during a typical binary's evolution because gravitational radia- 
tion is also produced during the hardening phase, as well as when the binary coalesces 
and rings down. 

Under this assumption, the maximum rest frame frequency f r , occurs for a circular 
orbit at 3 Schwarzschild Radii: 

, c 3 (M1+M2) 1 / 2 

Jrmx- UJkG W 

where Mi, and M2 are the black hole masses, G is the gravitational constant, and c is 
the speed of light. The minimum/,- is: 

_ 1 G(M l +M 2 ) 1 / 2 

J min — 3 \^) 

where a gw is the binary separation where gravitational radiation dominates, expressed 
by: 



^ = 0.0014 ^y ]■/«,,■/« (3) 



where t^ 1 / 4 is the coalescence timescale in Gyr. The change in the rest frame fre- 
quency is: 

f _ 96 f 5/3 G 3 MiM 2 (Mi +M 2 ) GjMi + M 2 ) 1 _ 1/3 
h ~ 5 h c 5 a 3 [ K 2 J { ) 

So, for each merger at a comoving distance d(z), gravitational radiation is emitted 
with a characteristic strain, h(f r 
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FIGURE 4. 4a-up left- shows the average observed characteristic gravitational strain amplitude as a 
function of observed redshift. Distance is measured from the center of the primary. Gravitational strain 
amplitude increases toward smaller redshifts as IMBHs sink into gravitational potential of the primary. 
This increase is rapid since the ratio of IMBHs binaries increases too. 4b-right- LISA sensitivity curve 
(black) and the total gravitational strain amplitude at different redshifts. The increase in mass of merging 
black holes with redshift results in increasing total gravitational strain amplitude and shifting toward 
higher observed frequencies. 4c and 4d on the bottom - gas accretion added. 
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